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1.  Introduction 

A correct  description  of  the  behavior  of  defects  in  thin-film  and  nanostructured 
systems  needs  to  take  into  account  both  the  strong  influence  of  outer  and  inner 
interphase  boundaries  and  the  fact  that  the  size  of  a defect  core,  where  classical 
solutions  are  singular  and  incorrect,  occurs  in  the  same  order  as  the  characteristic 
length  of  the  system  (film  thickness,  grain  size,  etc.)-  To  take  into  account  the 
influence  of  interphase  boundaries,  many  solutions  of  boundary  value  problems  for 
defects  have  been  obtained  in  the  framework  of  the  classical  theory  of  elasticity 
(see  [1-5]  for  a review).  To  solve  the  second  problem,  some  approaches  have 
been  proposed  which  have  been  aimed  at  dispensing  with  classical  singularity 
of  elastic  fields  within  defect  cores  (see  [6,  7]  for  a review).  The  present  paper 
represents  a brief  description  of  our  recent  results  [6-10]  dealing  with  nanoscale 
elastic  fields  within  and  near  cores  of  disclinations  [6,  7]  and  dislocations  [6-10] 
in  the  framework  of  gradient  elasticity.  The  main  result  shown  there  was  an 
elimination  of  displacement,  strain,  stress  and  energy  singularities  at  the  defect 
line.  It  is  worth  noting,  that  previous  continuum  models  for  such  kind  of  defects 
which  have  taken  into  account  couple  stresses  or  non-locality  (see  [6,  7]  for  a 
review),  do  not  dispense  with  the  singularity  in  the  displacement  or  strain  field, 
even  though  some  of  them  [11-13]  claim  elimination  of  stress  singularity. 

We  have  started  with  a simple  gradient  modification  of  the  linear  theory  of 
elasticity  in  the  form  [14] 

er  = A(tr  e)I  + 2pe  — cVz  [A(tr  e)I  + 2/ze] , (1) 

where  A and  fi  are  the  Lame  constants,  cr  and  e are  the  stress  and  strain  tensors,  I 
is  the  unit  tensor,  V2  denotes  Laplacian  and  c > 0 is  the  gradient  coefficient.  Using 
(1)  or  similar  theory,  the  authors  of  [14-20]  have  demonstrated  the  elimination  of 
classical  singularity  from  the  solution  for  the  strain  field  at  the  crack  tip. 
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Encouraged  by  these  results,  we  first  employed  [8,  9]  the  same  gradient  theory 
described  by  (1)  to  consider  dislocations.  In  particular,  four  dislocation  config- 
urations (i.e.  a screw  dislocation,  an  edge  dislocation,  and  the  dipoles  of  such 
dislocations)  have  been  considered.  It  has  been  shown  that  in  the  case  of  screw 
dislocation  [8],  the  elastic  strain  is  zero  at  the  dislocation  line  and  achieves  a maxi- 
mum value  (ss  ±bz/10Tr<Jc),  at  a distance  ~ \fc  from  it.  It  is  worth  noting  that  for 
an  atomic  lattice,  the  gradient  coefficient,  c,  can  be  estimated  [14]  as  \fc  « a/4, 
where  a is  the  lattice  constant.  With  the  Burgers  vector,  bz,  taken  to  be  equal  to 
a,  it  follows  that  the  aforementioned  maximum  value  is  estimated  as  « 12%.  In 
the  case  of  edge  dislocations  [9, 10],  we  have  also  shown  that  all  strain  components 
are  equal  to  zero  at  the  dislocation  line  achieving  maximum  values  (3-14)%  within 
the  dislocation  core  (r  < 4v/c). 

It  has  also  turned  out  for  both  types  of  dislocations  that  beyond  the  dislocation 
core  (ro  ~ 4,/c)  the  classical  and  gradient  solutions  coincide.  Furthermore,  it  has 
been  shown  that  the  total  displacement  varies  smoothly  across  the  dislocation  line 
in  contrast  to  the  classical  solution  which  suffers  a jump  there.  In  addition,  some 
of  the  displacement  components  which  are  singular  in  the  classical  theory  now 
become  finite  (edge  dislocation).  We  have  also  considered  the  displacement  and 
strain  distributions  for  dislocation  dipoles  and  shown  that  the  values  of  relative 
displacements  of  the  cut  surfaces  depend  on  the  dipole  arm,  in  contrast  to  the 
classical  theory  where  these  displacements  are  always  the  same.  As  a result,  two 
characteristic  distances  appear  naturally  in  this  approach:  ro  w 4 y/c  which  may 
be  viewed  as  the  radius  of  dislocation  core  and  do  » 10\/c  which  may  be  viewed 
as  the  radius  of  strong  short-range  nanoscale  interaction  between  dislocations. 

In  [6],  we  have  considered  the  elastic  fields  of  disclinations  within  the  gradient 
theory  described  by  (1).  We  have  found  the  elastic  strains  for  all  kinds  of  straight 
disclinations,  and  examined  the  interactions  between  disclinations  in  dipoles.  It 
has  been  shown  that  the  main  features  of  the  gradient  solutions  are  very  similar 
to  the  case  of  dislocations:  the  singularities  at  the  disclination  lines  are  eliminated 
from  the  strain  fields  which  are  equal  to  zero  or  attain  finite  values  there.  The  non- 
vanishing values  depend  strongly  on  the  dipole  arm,  d,  and  exhibit  a regular  and 
monotonous  (in  the  case  of  wedge  disclinations)  or  non-monotonous  (in  the  case 
of  twist  disclinations)  behavior  for  short-range  nanoscale  (d  < 10\/c)  interactions 
between  disclinations.  When  the  disclinations  annihilate  (d  — ► 0),  the  elastic 
strains  tend  to  zero  value.  Far  from  the  disclination  line  (r  10\/c)  gradient  and 
classical  solutions  coincide.  When  d \fc,  the  elastic  fields  of  a dipole  of  wedge 
disclinations  transform  into  the  elastic  fields  of  an  edge  dislocation  [9]  as  is  the 
case  in  the  classical  theory  of  elasticity. 

It  is  important  to  emphasize,  however,  that  in  the  framework  of  the  gradient 
elasticity  theory  described  by  (1),  the  stress  fields  of  dislocations  [8,  9]  and  discli- 
nations [6]  remain  as  in  the  classical  theory  of  elasticity  (same  as  in  the  case  of 
crack  problems  [14-17,  20]).  In  order  to  eliminate  the  singularities  also  from  the 
elastic  stresses  of  defects,  we  used  in  [7,  10]  a more  general  version  of  gradient 
elasticity  theory  which  has  been  also  utilized  by  Ru  and  Aifantis  [21]  (see  also 
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[22]).  The  constitutive  equation  of  this  theory  reads 

(l  - ci  V2)  a = (l  - c2V2)  [A(tr  e)7  + 2 fie] , (2) 

with  two  different  gradient  coefficients  ci  and  c2.  In  [21]  a rather  simple  mathe- 
matical procedure  analogous  to  the  one  contained  in  [15]  has  been  outlined  for  the 
solution  of  (2)  in  terms  of  solutions  of  classical  elasticity  for  the  same  boundary 
value  problem.  In  fact,  it  is  easily  established  (see  [15],  also  [6-10])  that  the  right 
hand  side  of  (2)  for  the  case  of  ci  = 0,  gives  the  classical  solution  for  the  stress  field 
which  we  denote  here  by  <r°,  while  the  solution  for  the  displacement  is  determined 
through  the  inhomogeneous  Helmholtz  equation  given  by 

(l  — c2V2)  u = u°  , (3) 

where  u°  denotes  the  solution  of  classical  elasticity  for  the  same  traction  boundary 
value  problem.  Eq.  (3)  implies  a similar  equation  for  strain,  e,  of  the  gradient 
theory 

(l  — c2V2)  e = e°  , (4) 

in  terms  of  the  strain  e°  of  the  classical  elasticity  theory  for  the  same  traction 
boundary  value  problem.  With  the  displacement  or  strain  field  thus  determined 
(which  is  obviously  independent  of  whether  c\  = 0 or  Ci  / 0),  it  follows  that  the 
stress  field  tr  of  (2)  can  be  determined  (for  the  case  ci  ^ 0)  from  the  equation 

(l-CiV2)cr  = £r°  , (5) 

where  cr°  denotes  the  solution  obtained  for  the  same  boundary  value  problem 
within  the  classical  theory  of  elasticity. 

Thus,  in  order  to  solve  equation  (2),  one  can  solve  separately  equations  (4) 
and  (5)  by  utilizing  the  classical  solutions  e°  and  cr°  provided  that  appropriate 
care  is  taken  for  the  extra  (due  to  the  higher  order  terms)  boundary  conditions  or 
conditions  at  infinity.  For  dislocations  and  disclinations,  this  problem’s  solutions 
are  accounted  for  by  assuming  that  the  strain  and  stress  fields  at  infinity  have 
the  same  characteristic  features  for  both  the  gradient  and  classical  theory.  The 
approach  has  already  been  applied  [21]  to  the  cases  of  screw  dislocations  and  mode- 
III  cracks  where  the  asymptotic  solutions  at  the  dislocation  fine  and  crack  tip  have 
been  found  demonstrating  the  elimination  of  both  strain  and  stress  singularities 
there. 

Below  we  report  the  exact  analytical  solutions  of  (2)  obtained  for  straight 
dislocations  and  disclinations. 


2.  Dislocations 

2.1.  CLASSICAL  SOLUTION 

Consider  a mixed  dislocation  whose  fine  coincides  with  the  z-axis  of  a Cartesian 
coordinate  system.  Let  its  Burgers  vector  be  b = bxex  + bzez  thus  determining 
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the  edge  (6X)  and  screw  (bz)  components.  In  the  framework  of  classical  elasticity 
theory,  the  total  displacement  field  is  described  by 

«°  = bxex^zBz  {arctan  | + |sign(y)[l  - sign(ar)]} 

+ Mr^{e-?-e»[(1'2,'>lnr+S]}’  (6) 

where  v is  the  Poisson  ratio,  r2  = x7+y2.  Here  we  use  a single- valued  discontinuous 
form  suggested  by  de  Wit  [23].  The  elastic  strain  field  £?•  reads  [23,  24]  (in  units 
of  l/[47r(l  - i/)])  by 

e°xx  = ~bxV  [(1  “ 2*')r2  + 2x2]/r4,  e°yy  = -bxy  [(1  - 2 v)r2  - 2 x2]/r4, 

£°Xy  = bxx{x2  - y2)/r4,  e°xz  = -bz(  1 - %/r2,  e°yz  = bz(  1 - u)x/r\  (7) 

and  elastic  stress  field  cr?-  is  [23,  24]  (in  units  of  /z/[2tt(1  — i/)]) 

= £°xx(V  = °)»  *lv  = £°vv(.V  = °)>  a°zz  = "("xx  + 0> 


<r°  =e° 

uxy  °xyi 


vv  ^yy^  ”zz  — ^v'xx  ^ "yyh 

°lz  = e°xz>  =e°  (8) 


Eyzi 


Fields  (6)  (y-component),  (7)  and  (8)  are  singular  at  the  dislocation  line. 

The  elastic  energy  W°  of  the  dislocation  per  unit  dislocation  length  is  [24] 


W°  = T-  (bl  + -^-)\n- , 
4tt  V 1 - v J r0 


where  R denotes  the  size  of  the  solid  and  ro  is  a cut-off  radius  for  the  dislocation 
elastic  field  near  the  dislocation  line.  When  ro  ->  0,  W°  becomes  singular. 


2.2.  GRADIENT  SOLUTION 


Let  us  now  consider  the  corresponding  dislocation  fields  within  the  theory  of  gra- 
dient elasticity  given  by  (2).  As  described  in  Section  1,  one  can  obtain  the  solution 
of  (2)  by  solving  separately  equations  (3)-(5).  They  can  be  solved  [9]  by  using  the 
Fourier  transform  method.  Omitting  intermediate  calculations,  we  give  here  only 
the  final  results.  For  the  total  displacements,  solution  of  (3)  gives  [8,  9] 

u = V°  ~ 47r(i-  t/)  {te*  2xy  + ev  fa2  ~ 7-2 + e»$°} 

+ sign (y)  r ip4e~M^ds  , (10) 

2tt  J0  J.  + S 2 

where  u°  is  given  by  (6),  4>o  = (l-2i/) K0(r/y/c^),  4>2  = [2 c2/r2 -K2{r / y/c^/r4, , 
Kn{r / <yc2)  is  the  modified  Bessel  function  of  the  second  kind  and  n = 0,1,... 
denotes  the  order  of  this  function. 
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For  the  elastic  strain,  solution  of  (4)  gives  [8,  9]  £y  = + efj,  where  £?•  are 

given  by  (7)  and  efj  (in  units  of  l/[27r(l  - i/)])  by 

«£=  b*y  l(y2  - + (3x2 *  - 2/2 *)4>2],  eff  = bz  (1  - v)yr2$i/2  , 

= 6*2/  [(®2  “ vr2)$i  - (3a:2  - y2)§2],  e9yrz  = - bz  (1  - u)xr2§ i/2  , 

= -&**  b2$x  + (a:2  - 3y2)$2],  (11) 

where  4>i  = Ki  (r/A/c2)/(N/c2r3).  For  the  stresses,  the  solution  of  (5)  gives  [10] 
<ry  = 0y  + ofj,  where  Oy  are  given  by  (8)  and  off  (in  units  of  ju/[ar(l  — i/)])  by 

<x  = = 0, c2  cx),  = e£(«/  = 0,c2o  Cl),  off  = o(off  + off), 

off  = exp(C2  **  Cx),  Off  =£®f(c2  •<->•  Cx),  Off  = £ff(c2  -HCx).  (12) 

The  main  feature  of  the  solution  given  by  (10)-(12)  is  the  absence  of  any 
singularities  in  the  displacement,  strain  and  stress  fields.  In  fact,  when  r 0,  we 
have  K0(r/y/cj;)\r^0  -»  -7  + In (2 y/cj^/r),  K^r/y/c^)  -»  y/c^/r,  K2(r/y/c£)  -> 
2 Cfc/r2  — 1/2,  (A:  = 1,2)  and,  thus,  uy  is  finite,  £y  ->  0,  oy  — ► 0.  The  fields  of 
displacements  (10)  and  strains  (11)  have  been  discussed  in  detail  in  [8,  9]  within 
a special  version  of  gradient  elasticity  theory  (cx  = 0).  The  stress  fields  (12)  have 
been  examined  in  [10].  It  has  been  shown  there  that  they  attain  their  extreme 
values  [\oxx\  « 0.45 /r,  |oTO|  ns  \axy\  ns  0.27 fj,  and  |ay|  ns  0.25 for  bi  = a = 4 y/ci[ 
and  v = 0.3)  at  a distance  r:  a/4  from  the  dislocation  line.  Far  away  from  the 
dislocation  core  (r>r0  ns  4 y/ck),  the  gradient  solutions  coincide  with  the  classical 
ones  [8-10]. 

Using  (12),  one  can  find  the  elastic  energy  of  the  dislocation  within  the  gradient 
elasticity  given  by  (2)  as  follows  [10] 

<13> 

where  7 = 0.57721566 ...  is  Euler’s  constant.  Thus,  we  obtain  a strain  energy 
expression  which  is  not  singular  at  the  dislocation  line. 


3.  Disclinations 

3.1.  CLASSICAL  SOLUTION 

Consider  a disclination  of  general  type  with  Frank  vector  us  = toxex  +u)yey  +u)zez  in 

an  infinite  elastic  medium.  The  scalars  ux  and  u>y  determine  the  twist  components 

of  the  disclination  while  uz  determines  its  wedge  component.  Let  its  line  coincides 

with  the  z-axis  of  the  above  coordinate  system.  For  such  an  isolated  disclination, 

both  classical  and  gradient  solutions  themselves  have  no  physical  meaning  because 

they  are  not  screened  but  they  may  be  used  in  modeling  screened  disclination 
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configurations  as  basic  elements.  The  classical  solution  for  elastic  strain  fields  sp- 
reads [23]  (in  units  of  1/[47t(1  — i/)]) 

2vr2  - x2  +y2  2vr2-Zx2-y2  |\  y2\ 

■ —4 +UJvzV -4 + W,|(l-2v)lnr+^j, 


'VV 


-xy 


2vr 2 - x2  - 3 y2  2vr2  + x2  - y2 

U)xzx 3 +U)szy 2 — 


+ uz  |(l-2i/) lnr  + ^ j , 
= (“vx  ~ uxy)zX  ~ri  V - wz^,  e°xz  =wy~  —wx  |(1  — 2i/)lnr  + ^ | , 


£yz  — r2 


— wy|(l  — 2i/)lnr+  J , 


(14) 


and  for  the  stress  fields  it  may  be  written  (in  units  of  h/[2-k(1  — i/)])  as 


r°  — 1 1 m „o  _ _0  _ n\  „o  _ _o  _o  _ _o  _ „o 

£xx\v  ^)i  ®yy  £yy\V  °zt/  £xi />  uxz  ~ £xz>  ayz  £yzi 


'XX  ^xx\ 


xy  ^xy> 


yz  i/z’ 


n •i'  y 

(jzz  = -uix  z 2v-^  — u)y  z 2v-^  + wz  2v  In  r . 


° pi. 

Most  of  the  components  in  (14)  and  (15)  contain  singular  terms  ~ In  r. 


(15) 


3.2.  GRADIENT  SOLUTION 


The  gradient  solutions  have  been  originally  obtained  for  a disclination  dipole 
within  both  the  gradient  theories  considered  in  Section  1.  Solving  (4),  we  solve 
finally  for  an  individual  disclination  under  consideration  the  strain  field  [6]  £y  = 
Sy  + efj,  where  e^-  are  given  by  (14)  and  efj  (in  units  of  l/[47r(l  — i/)])  by 


egr 

°xx 

= 

u}x  2 xz  {(y2  - vr2)$i  + (x2  - 

■3  j/2)$2} 

+ 

oJy  2yz  {{y2  - vr2)$i  + (3x2 

- j/2)$ 2}  + ujz  {$0  + r2{x2  - 

-2/2)*2}. 

egr 

°yy 

= 

ux  2 xz  {(x2  - vr2)$i  - (x2  - 

-3  y2)$2} 

+ 

uiy  2 yz{(x2  - z/r2)$i  - (3r2 

- 2/2)*2}  + uz  {$0  - r2(ar2  - 

-y2)$2}, 

egr 

xy 

= 

ux2yz{-x2$1  + (3r2  -y2)$2} 

- 

Uy  2 xz  {y2$i  + {x2  - 3t/2)$2 

}+w2  2xyr2$2, 

egr 

°xz 

= 

Wi  {-$o  - r2(a:2  - y2)$2}  - 

u)y  2xyr2§2, 

egr 

cyz 

= 

u)y{-$o+r2{x2  -y2)$2}  - 

ux  2xyr2$2, 

(16) 

where 

are  the  same  as  in  Section  2 

.2.  For  the  stress  field,  the  solution  of 

(5)  gives  [7]  Oij  = cr?-  + where  <7?-  are  given  by  (15)  and  <r?r  (in  units  of 

MM1  - ")])  by 

axx  = £xrx(l/  = 0,c2<+  Cl),  a9yry  = egry{ V = 0 ,c2  ci), 

<t9/z  =2v{(ajxx  + (jyy)zr2$1(c2  Ci)  + wz$0(v  = 0,c2  «->ci)}, 
^y=£xry(c 2-W-Ci),  0*TZ  = egrz(c2  O Ci),  <=^(c2«Ci).  (17) 
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Using  the  limiting  transitions  noted  in  Section  2.2,  it  is  easy  to  show  the 
total  elimination  of  classical  logarithmic  singularity  from  elastic  fields  (16)  and 
(17).  In  [6,  7],  we  have  considered  similar  elastic  fields  of  disclination  dipoles 
in  detail  and  found  that  they  are  equal  to  zero  or  attain  finite  values  at  the 
disclination  lines.  The  non-vanishing  values  depend  strongly  on  the  dipole  arm,  d, 
and  exhibit  a regular  and  monotonous  (wedge  disclinations)  or  non-monotonous 
(twist  disclinations)  behavior  for  short-range  nanoscale  ( d < 10-^/c)  interactions 
between  disclinations.  When  the  disclinations  annihilate  (d  ->  0),  the  elastic 
strains  and  stresses  tend  to  zero.  Far  from  the  disclination  line  (r  » 10a/c) 
gradient  and  classical  solutions  coincide.  When  d -C  \fc,  the  elastic  fields  of  a 
dipole  of  wedge  disclinations  transform  into  the  elastic  fields  of  an  edge  dislocation 
[9,  10]  as  is  the  case  in  the  classical  theory  of  elasticity. 


4.  Conclusions 

The  gradient  elasticity  described  by  (2)  has  been  employed  in  the  consideration  of 
nanoscale  short-range  elastic  fields  and  interactions  of  dislocations  and  disclina- 
tions. Exact  analytical  solutions  for  the  displacements,  strain  and  stress  fields  and 
elastic  energies  of  dislocations  have  been  reported  which  demonstrate  the  elim- 
ination of  any  singularity  from  the  elastic  fields  and  energies  at  the  dislocation 
line.  For  disclinations,  the  gradient  solutions  have  been  given  for  the  strain  and 
stress  fields  where  the  classical  singularities  at  the  disclination  lines  also  disap- 
pear. These  new  non-singular  elastic  fields  are  considered  as  especially  useful  for 
modeling  nanoscale  behavior  and  interactions  of  dislocations  and  disclinations  in 
thin-film  or  nanostructured  solids. 
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